Kondo Tunneling into a Quantum Spin Hall Insulator 
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The physics of a junction composed of a normal metal, quantum dot and 2D topological insulator 
(in a quantum spin Hall state) is elucidated. It maifests a subtle combination of Kondo correlations 
and quantum spin Hall edge states moving on the opposite sides of the 2D topological insulator. 
In a narrow strip geometry these edge states interact and a gap opens in the edge state spectrum. 
Consequently, Kondo screening is less effective and that affects electron transport through the junc- 
tion. Specifically, when edge state coupling is strong enough, the tunneling differential conductance 
develops a dip at zero temperature instead of the standard zero bias Kondo peak. 

PACS numbers: 71.10.Pm,73.43.-f,72.15.Qm,73.23.-b 



Background: Topological insulators (TI) and topolog- 
ical superconductors (TS) have recently become a sub- 
ject of intense theoretical and experimental studies.— ~— 
In most of these studies, properties of topological states 
are analyzed using classification of the electronic phases 
according to the pertinent topological invariants r^— 

Two-dimensional TI arc systems in which time reversal 
symmetry is respected but spin rotation invariance is vio- 
lated. They differ from ordinary insulators as they posses 
a non-trivial topological number Z-2 implying a gapped 
spectrum in the bulk and gapless spectrum of quantum 
spin Hall (QSH) edge states^ (or helical modes). A 
helical mode consists of a Kramers degenerate pair of 
states propagating in opposite directions along the same 
edge. Elucidating the QSH helical modes in experiments 
is rather elusive since both edge states moving on op- 
posite boundaries contribute to the conductance^ This 
obstacle can be alleviated if there are strong correlations 
of edge states on the opposite edges, leading to the emer- 
gence of massive Dirac fcrmions instead of the massless 
ones that prevail in the absence of such correlations. 
Motivation: In order to employ this feature, we note 
that the physics of QSH state can be studied in tun- 
neling experiments In particular, we focus on tun- 
neling through an interacting quantum dot tuned to the 
Kondo regime. In the absence of correlations between 
edge states, Kondo tunneling into a TI with gapless he- 
lical edge state spectrum displays a finite zero bias dif- 
ferential conductance at zero temperature^, as is the 
case for tunneling into a metal. Once the coupling be- 
tween the helical states on the opposite sides of the TI is 
present, a gap opens in the edge states spectrum, which 
naturally affects the tunneling conductance^ Screening 
of the magnetic impurity becomes less effective and the 
conductance develops a dip at zero temperature. 
The main objective: The main goal of the present 
work is to substantiate the above qualitative analysis. It 
is argued that interaction between QSH edge states on a 
2D TI has a clear signature that turn it to be experimen- 
tally detectable as it leads to unusual behavior of Kondo 
tunneling at zero temperature. 

Main results: We have theoretically analyzed tunnel- 
ing through a device consisting of a tunnel junction com- 



posed of normal metal (NM) interacting quantum dot 
(QD) and 2D TI as shown in Fig. [Q The TI is modeled 
as a stripe with topological states moving on its oppo- 
site sides with a coupling between them. Interrelation 
between the Kondo physics prevailing in the QD and the 
QSH physics prevailing in the TI is then studied. Kondo 
tunneling is analyzed in the weak (T 3> Tk) and strong 
(T < Tk) coupling regimes, where Tk is a Kondo tem- 
perature. It is shown that in both regimes, the conduc- 
tance develops a dip when the temperature T decreases 
and crosses an energy threshold v specifying the coupling 
strength between the QSH edge states. Since this setup 
is experimentally feasible we argue that coupling between 
QSH edge states residing on opposite edges of a 2D TI 
has a definite experimental signature. 



FIG. 1: NM-QD-TI junction. Helical states exist on the edges 
of the TI. The QD can be gated to adjust the level position 
in the dot and/or the tunneling rates ti, and tn. 

Anderson Model for NM-QD-TI Junction: The 

junction consists of the following ingredients (see Fig. [TJ: 
the NM lead, held at bias voltage V; the TI with coupling 
between helical modes on its two opposite sides; the QD 
in a tunnel contact with the NM and with the TI. The 
Hamiltonian of the junction is then written as, 

H = H L +H R + H D +H T . (1) 

The first term on the right hand side of Eq. (p} is the 
Hamiltonian of the NM lead, 



Hr = 



E 
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The second term, Hr, is the low energy Hamiltonian of 
the 2D TI, 

v ar k a\ ar aknr + [ a lai a k<y2 + b..c.j,(3) 
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where v is the coupling strength between the helical 
fermions on the opposite sides (r = 1,2) of the TI, and 
v ar = (— l) r_1 (7u, with v being the Fermi velocity. 

The coupling v represents the effective mass of Dirac 
fermions. It is convenient to work in a basis for which 
the TI Hamiltonian ([3]) is diagonal, obtained through the 
canonical transformations, 



Oikal = Ukjkcr+ ~ PVkJkcr- 
aka2 = dVk r Jka+ + Uklko- 



(4) 



where 



£k : 



Ek + hvk 



2et 



Ek — hvk 



2e k 



^(hvk) 2 + v 2 . 
The Hamiltonian ([3]) then takes the form, 



H r , 



E 

ka,p— zb 



(5) 



The third term on the right hand side of Eq. ([T]), Hp, is 
the QD Hamiltonian, 



Hn = 



(6) 



Here d\ , d a are creation and annihilation operators for 
the dot electrons; n a = d' a d a , and U > is the strength 
of the Coulomb repulsion. The last term on the right 
hand side of Eq. ([1]), H T = Htl _|_ Htr, is composed 
of Htl, describing electron tunneling between the dot 
and the normal metal lead, and Htr describing electron 
tunneling between the dot and the TI. 

Denoting by ^x(r) and iPr(t) the electron field op- 
erators in the NM and the TI and assuming the dot is 
positioned at r — 0, we have, 



#TL=^^[vL(0K+h.c. 
Htr = E [^L(°K + h - c - 



(7) 



where t^,R are the corresponding tunneling amplitudes 
assumed to be spin independent. 

The Spin Hamiltonian: The quantum dot is tuned by 
gate voltage to have a single electron in its ground state. 
The depth |eo| of the dot level is assumed to be much 
larger than the tunneling width, 



r. 



T L = 27rf! PL , = 2nt 2 R p R , (8) 



where pl is the electron density of states in the NM, while 
pn is the electron density of states belonging to the gap- 
less part of the TI. Hence, charge fluctuations can be in- 
tegrated out using the Schricffcr- Wolff transformation^. 
It projects out zero and two electron states in the dot 
and maps the Hamiltonian H, Eq. ([T]) onto an effective 



spin Hamiltonian H acting in a subspace of states where 
there is one and only one electron on the dot. The effec- 
tive Hamiltonian in this subspace reads, 

h k = ^<*MKi> Vw(0), (9) 



where a, a' = L, R, and the matrix h^a' is defined as, 

^aa' = ~^ H aa 'S a{ j' -\- — J aa rS • T aa i . 

Here r is the vector of the Pauli matrices and S is the 
spin operator of the electron residing on the QD. The 
coupling constants are 

It t i It t i It t i It t i 
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Scaling Equations Using the poor man's scaling tech- 
nique, we obtain scaling equations for the dimensionless 
coupling constants j aa > = J aa i w 'p a p a ' ', 



djhh 
d\nD 

djRR 

dh\D 

djLR 



Ull + Jlr) 



— Urr + Jlr) 



= -JLR [JLL +JRR 



d\nD 

subject to the initial conditions, 

vCT^ / l i 



(10) 

(ii) 

(12) 



jaa> (D) = 



* V kol U-eo 



If v is the smallest energy scale, the solutions of these 
equations are, 



jaa'( T ) = 



%/TVlq^ 1 



r In (T/T K ) ' 
and the scaling invariant, (the Kondo temperature), is 



(13) 



T, 



K 



D exp < — 



r 



(14) 



When v is not the smallest energy scale, the scaling 
behavior is more complicated. The flow diagram still has 
a fixed point at infinity, but the Kondo temperature turns 
out to be a sharp function of i/ i 25 i 26 For v ^> Tjc, the 
scaling of j aa > (T) depends on whether the temperature 
is higher or lower than v. For T 3> v Tk, the gap in 
the edge state spectrum can be neglected and the scaling 
of jaa' (T) are given by Eq. (|13|) with the scaling invariant 
T K , Eq. ([H]). For v > T K > T the scaling of j LR and 
3rr terminates at D ~ v and for D < is we have just one 
scaling equation, leading to a fixed point at ]ll — > oo 
and a Kondo temperature that depends on v. However, 
we will see that the tunneling conductance diminishes in 
the low temperature regime, T <C v. Therefore in the 
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following discussions we may speak of a single Kondo 
temper asture, given by Eq. (| 14[) . 

Calculations of the Tunneling Conductance in the 
Weak Coupling Regime are carried out below using 
the Keldysh technique in order to treat a system out of 
equilibrium. The required quantities to be used elsewhere 
below are the Keldysh Green's functions (GF) 



- 9 L a sr 

g£ 



(15) 



where the subscript a refers to the normal metal (L) 
and TI (R) or the dot (f) electron GF while the super- 
scripts refer to retarded (R), advanced (A) and Keldysh 
(K) types of the GF. For a = L,R the explicit expres- 
sions are, 



Sl = ~9l 

-R,A 



iTT PL , g^(e) = -2inp L {l-2f(e)), 



9r (e) = -KRR x(e±if7), 
g* (e) = -2mp R Im X (e) (l-2/(e)), 

where /(e) is the Fermi function and 



X{z) 



(16) 



y/ Z 2 — v 2 



(17) 



The current operator for tunnelling from the NM to 
the TI is 

1 = |E{^L(0)^^(0)-h.c.}, (18) 



and the differential dc conductance is given by, 
where the expectation value of the current is 



T{V) = (Uk~ 1 T{Uk 1(0))). 



(19) 



In this expression, T is the time-ordering operator and 
Uk is the evolution operator under Hk , 



U K = Tcxp^ -- / dtH K {t)\. (20) 



To lowest (second) non- vanishing order of perturbation 
theory in the dimensionless parameters j aa i and k aa i = 
K aa i y/ p a p a ' , the conductance is 

2 

G 2 = {kl R + 3j 2 LR }w(T,V), (21) 
where W(T, V) = w(T, V) +w{T,-V), 

OO 

w{T, V) 



4T J cosh 2 



The third order correction to the conductance is 
37re 2 

a=L.R 



G, 



2h 



f LR W(T,V) J2 JctMD), 



(22) 



where D is a high energy cut off, 



£l{d) 



1 



In 



D 



\y/(v + eV) 2 +T 2 . 
D 



In 



Cr{D) 



In 



yj{v- eV) 2 + T 2 , 
D 



s/v 2 + T 2 



The third order correction G3 contains a logarithmic 
term which is, strictly speaking, not small. As a result, 
G3 is not small compared with G2 and the expansion up 
to terms cubic in j aa ' is insufficient. Instead, we derive 
the conductance in the leading logarithmic approxima- 
tion using the RG equations flTUJ) ™ (fT2")h 

In the following analysis we split the second order con- 
tribution to the conductance, Eq. (|2ip . in two parts: 
one part is due to the exchange cotunncling, which is 
proportional to j\ R \ and one part is due to regular co- 
tunneling, which is proportional to k\ R . G3, Eq. (|22[) . 
increases when the temperature and voltage are lowered, 
demonstrating the Kondo anomaly. The regular cotun- 
ncling contribution contains k\ R which does not grow at 
low temperatures and bias, and therefore it does not con- 
tribute to the Kondo effect. The exchange cotunneling 
contains term j\ R which demonstrates logarithmic en- 
hancement of the conductance at low temperatures [see 
Eq. (p~3|) ] and contributes to the Kondo effect. Therefore, 
we single out the exchange contribution in the second or- 
der term, 



G° xch (L>) 



37re 2 
Ah 



f LR (D)W(T,V). 



(23) 



The resulting condition of invariance of the conductance 
under the transformation, which corresponds to "poor 
man's scaling" , has the form, 



d 



d\nD 
3ire 2 



+ - 



21i 



G e 2 xch (D) 



f LR W{T,V) 



E 

a=L,R. 



j a aC a (D) \ - 0. (24) 



Within the accuracy of this equation, when differenti- 
ating the second term, we should neglect any implicit 
dependence on D through the couplings j a a'- Eq. d2U) 
yields the scaling equation (fT2l . 

The renormalization procedure should proceed un- 
til the bandwidth D is reduced to a quantity 
d(T, V,v) =Mx[T,eV>], At this point, the third order 
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correction to the conductance is much smaller than GSf 1011 
and the current and conductance can be calculated in 
the Born approximation, as in Eq. (|23|) . This situation 
is similar to the problem considered in Ref^i, but the 
coupling between the helical modes in the TI introduces 
an additional energy scale v which affects the Kondo 
physics. When D reduces below v, renormalization of 
both 2lr and ]rr stops. Then, taking 



d(T,V,v) = ^u 2 + {eV) 2 +T 2 , 
we get an expression for the conductance for T, |eV|>Xx, 
3n 2 G W(T, V) 



G(T, V) 



where 



Go 



16 In 2 (d{T,V,p)/T K ) 



(25) 




T/T K 



FIG. 2: Conductance as function of temperature in the weak 
coupling regime. The curves a, b, c and d correspond to 
v = 0, 5Tk, WTk and 15Tk, respectively. The case v = 
corresponds to the standard Kondo effect in NM-QD-NM 
junction.— 
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FIG. 3: The nonlinear conductance (I25[) as a function of ap- 
plied bias for T = 5Tk- The curves a, b, c and d correspond 
to v = 0, 5Tk, WTk and 153V, respectively 

The total differential conductance is displayed in Fig. 
d]for V = 0. It is seen that when the gap in the spectrum 



of the TI exceeds Tk (see curves b, c and d), the conduc- 
tance increases when the temperature is lowered, reaches 
its peak at finite temperature and begins to decrease with 
further lowering the temperature. This occurs because 
electron needs an extra energy in order to covercome the 
gap v between the Fermi level of the left lead and the bot- 
tom of the conduction band of the TI. This energy can be 
transferred from the internal energy of the electron liq- 
uid in the lead and/or from the voltage difference across 
the barrier. When the temperature is lower than the 
gap, the zero bias conductance decreases with temprera- 
ture. For comparison, it is noticed that the conductance 
for the gapless spectrum of the TI (curve a) increases 
when the temperature is lowed which corresponds to the 
standard feature of Kondo tunneling through the NM- 
QD-NM junctional 

The coupling between the helical states manifests it- 
self also in the nonlinear tunneling conductance. Fig. [3] 
displays the nonlinear conductance as a function of the 
applied voltage for T = 5Tk and for several values of 
v. The zero bias peak in the conductance correspond- 
ing to the standard Kondo tunneling decreases with the 
coupling energy v (curves a and b) and splits into two dis- 
tinct peaks when v exceeds the temperature T (curves c 
and d) . This is one of the central results of the present 
paper since it combines the Kondo and the QSH physics. 
Conductance Calculations in the Strong Coupling 
Regime: When T < Tk, the mean field slave boson ap- 
proximation (MFSBA) is used to calculate the zero bias 
tunneling conductance. In the limit U — > oo, the dot can 
be empty or singly occupied. The dot operators are writ- 
ten as d a = f a and d a — f a b where the slave fermion 
operators f a and the slave boson operator b satisfy con- 
straint condition, 

Q = Y,ftf- + ^b=l. 

a 

This condition is encoded by including a Lagrange mul- 
tiplier A in the total action S. In the mean field approx- 
imation, we replace the Bose operators b and by their 
expectation values, 



b^b Q , tf^b , 



b 



(bfb). 



At the mean field level the constraint condition is satis- 
fied only on the average. 

The partition function Z(a) is calculated by integrat- 
ing the action over slave fermion field. Here the source 
field a is coupled to the current operator, 



I = 



iet L 

h 



-h.c. 



(26) 



The effective action in the MFSBA is Gaussian and de- 
pends on two real numbers, the boson field bo and the 
chemical potential (Lagrange multiplier) A. Integrating 
the action leads to the partition function, 



\nZ(a) = -tr hi jcj 1 



ea q t 2 L b 2 



[g L ,r x ]}, (27) 
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where 



GJ L = gj L - tib^L - t R b z g R . 



(28) 



Here gf is the GF of the (non-interacting) electron in the 
QD with shifted energy level, eo — > £q = £ o + A, 



9 f /A (e) = 



1 



e - e ± 



2^(1-2/(6)) 
(e-e ) 2 + ?7 2 ; 



where Tjf is the Kondo temperature given by Eq. (|14l) 
Then the linear conductance for T < is 



G = 



ttT k Go f deS f (e) {T L + T R )Im X (e) 



AT 



cosh 2(_L) T L +T R Im X (e) 



(31) 



where 




FIG. 4: The zero bias conductance as a function of temper- 
ature in the strong coupling limit (T < Tk)- The curves a, 
b, c, d and e correspond to v = 0, 0.2Tk, OATk, Q.GTk and 
0.8Tk, respectively. Here we take Tl = Fr. 

The MFSBA is reliable in equilibrium, V = 0. There- 
fore we will consider below the temperature dependence 
of the zero bias conductance. In equilibrium, the mean 
field solutions for bo and A minimize the free energy, 

F = -Tj2 tT lnG7 1 (iw„) + Xb 2 , 



where the last term is the slave boson kinetic part of 
the free energy due to the constraint, Gj (iui n ) is the 
Matsubara's GF. The mean field equations are 



&o = — arctan ( - [Vj^ — - 
7r V b o l 



A 



ttD 

2 




These equations are solved for bo and A with the solu- 
tions, 



ttT, 



K 



£o + A = 



7r 3 (T 



2T 



(30) 



Sf{e) 



7(e) 



TKl{e) 



{e-loY + {T K Y^{e)' 
7T r L + T R Imx(e) 
2 



Go is given by Eq. (|26l) . x(e) is given by Eq. (|T7j) . 

The zero bias conductance as a function of temperature 
is shown in Fig. Hfor v = 0, 0.2T K , 0AT K , 0.6T K and 
0.8Tk- The temperature dependence of the conductance 
is similar to the one shown in Fig. [5] for T > Tk- In both 
cases, when T > is, a peak in the conductance occurs 
(see curves b, c, d and e). This peak decreases with v 
and shifts towards high temperatures. For T < v, the 
conductance decreases with decreasing in T and vanishes 
when T — > 0. This is a manifestation of the coupling 
between the QSH states on the opposite sides of the TI. 
For v = 0, the conductance has a usual peak at T = 0. 

Summary: The linear and non-linear conductance in a 
system consisting of a quantum dot (tuned to the Kondo 
regime) connected to a metal lead on the left side and 
to a topological insulator on the right side is evaluated 
using Kcldysh technique. Rcnormalization group analy- 
sis is performed in the weak coupling regime (T ^> Tk) 
while the MFSBA is used at at the strong coupling regime 
T < Tk- When the coupling energy v between QSH he- 
lical states exceeds the temperature, the differential con- 
ductance develops a dip at low temperature, in contrast 
with the standard zero bias anomaly Kondo peak prevail- 
ing in QD connected to metallic leads. Our analysis then 
shows that the Kondo resonance is very sensitive to edge 
state coupling and can be used to probe the features of 
QSH helical modes. 
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